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Abstract 

We present the ZEN (z equals nine) survey: a deep, narrow J-band search for 
proto-galactic Lya emission at redshifts z ~ 9 . In the first phase of the survey, 
dubbed ZEN1, we combine an exceptionally deep image of the Hubble Deep Field 
South, obtained using a narrow band filter centred on the wavelength 1.187/im, with 
existing deep, broad band images covering optical to near infrared wavelengths. 
Candidate z ~ 9 Lya-emitting galaxies display a significant narrow band excess 
relative to the J s -band that are undetected at optical wavelengths. We detect no 
sources consistent with this criterion to the 90% point source flux limit of the NB 
image, Fnb = 3.28 x 10 -18 ergss -1 cm -2 . The survey selection function indicates 
that we have sampled a volume of approximately 340 /i~ 3 Mpc 3 to a Lya emission 
luminosity of 10 43 /i -2 ergss -1 . When compared to the predicted properties of z ~ 9 
galaxies based upon no evolution of observed z ~ 6 Lya-emitting galaxies, the 'vol- 
ume shortfall' of the current survey, i.e. the volume required to detect this putative 
population, is a factor of at least 8 to 10. We also discuss continuing narrow J-band 
imaging surveys that will reduce the volume shortfall factor to the point where the 
no-evolution prediction from z ~ 6 is probed in a meaningful manner. 



1 Introduction 



Determining during which epoch and under what conditions the first galaxies 
and stars formed is a key goal of galaxy evolution studies. Observations of 
distant galaxies currently extend to redshifts z ~ 7 (Kneib et al. 2004). In 
addition, a number of studies have demonstrated that selecting faint galaxies 
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displaying i— z colours characteristic of the Lyman Break signature at redshifts 
z ~ 6 is an effective method to identify significant (> fO galaxies) samples of 
high-redshift galaxies (Bouwens et al. 2004a; Dickinson et al. 2004; Stanway 
et al. 2004; Malhotra et al. 2005). The addition of a narrow band (NB) filter 
to such broad band observations robustly identifies Lya-emitting galaxies at 
redshifts that place the emission line in the filter bandpass (Hu et al. 2004; 
Kodaira et al. 2003; Rhoads et al. 2004). In the case of Lya-emitting galaxies 
at z > 5, NB filters can exploit narrow spectral regions free from bright night 
sky emission lines. The main benefit of this approach is that the spectral 
location and profile of the Lja emission line can be identified in a relatively 
straightforward manner - leading to samples of galaxies with very high spectral 
completeness (e.g. Hu et al. 2004). The main drawback of course is that the 
narrow spectral region viewed by the NB filter greatly reduces the volume 
sampled compared to a broad band survey of comparable depth and field. 

The properties of individual z ~ 6 galaxies now provide relatively detailed 
constraints upon the the very earliest epochs of star formation (SF) - point- 
ing to an earlier epoch of very intense SF just beyond the redshift limit of 
current, optical surveys. Two of the brightest galaxies known at z ~ 6 have 
been detected at 3.6yum and 4.5/im using the Spitzer infrared space telescope 
(Eyles et al. 2005). Analysis of the rest-frame UV to optical spectral energy 
distribution (SED) of these galaxies indicates the presence of stellar popula- 
tions of masses ~3 x 10 10 M and ages ~ 400 Myr. Within the cosmological 
model assumed in these proceedings (see below) the epoch corresponding to 
the onset of SF in these galaxies is z ~ 9. 

A key factor that determines the visibility of Lya-emitting galaxies at z > 6 is 
the ionisation state of the inter-galactic medium (IGM). As the neutral fraction 
of uniformly distributed hydrogen gas in the Universe increases beyond 1 part 
in 10 5 the associated optical depth at rest wavelengths blueward of 1216A 
exceeds unity (Gunn and Peterson 1965). Observations of "dark regions" - 
- consistent with an IGM optical depth r > 1 - in redshift z > 6.2 QSOs 
suggest that we may be witnessing the onset of neutrality in the IGM at 
these redshifts (Becker et al. 2001; Fan et al. 2002). As the IGM neutral 
fraction increases further, the absorption profile develops a strong damping 
wing extending to rest wavelengths A > 1216A - potentially absorbing the 
Lya feature in high-redshift SF galaxies. However, the observation of a Lya- 
emitting galaxy at redshift z = 6.56 (Hu et al. 2002) is not inconsistent with 
observation of a possible GP effect in lower redshift QSOs. Haiman (2002), 
Santos (2004) and Barton et al. (2004) describe the ionising effect of a star 
forming galaxy embedded in a neutral IGM and note that, depending upon 
the exact assumptions made regarding the mass and star formation properties 
of the source and the physical conditions present in the IGM, the galaxy will 
form a local HII region of sufficient size to permit transmission of a partially 
attenuated Lya line and associated continuum. 
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The visibility of Lya emission in high-redshift SF galaxies is important as, 
although z > 7 galaxies will appear as continuum drop-out sources, e.g. z — J 
(Bouwens et al. 2004b), confirming the redshift of such faint candidates (H > 
27) on the basis of continuum features observed in deep near infrared (NIR) 
spectra is challenging to the point of being impractical. Therefore, a NIR NB 
survey for Lya-emitting galaxies at z > 7 has the potential to detect Lya- 
bright galaxies for which obtaining a spectroscopic redshift via follow-up NIR 
spectroscopy using 8-10m class telescopes is a realistic goal. In the remainder 
of these proceedings we describe a dedicated search for high-redshift star 
forming galaxies, employing an extension of broad and narrow-band selection 
techniques applied at optical wavelengths to the NIR wavelength regime. In 
particular we focus upon the application of a narrow J-band filter centred at 
A = 1.187/im to detect the signature of Lya emitting galaxies located about 
a redshift z = 8.8 (termed z ~ 9 in the following text). 

In these proceedings we adopt a cosmological model described by the param- 
eters Q M = 0.3, Q A = 0.7, h = Ho/100 kms^ 1 Mpc" 1 = 0.7. All magnitudes 
are quoted using the AB system. 



2 Constructing the experiment: going deep and rejecting interlop- 
ers 



We assume that a redshift z ~ 9 Lya-emitting galaxy will display a signifi- 
cant narrow band excess relative to the J-band, in addition to displaying a 
continuum break consistent with almost complete attenuation of photons at 
rest frame A < 1216A. In order to generate an effective survey for such sources 
certain additional factors must be considered: 

(1) NIR continuum imaging data must achieve a limiting depth of AB > 25.5. 
The brightest emission line galaxies confirmed at redshifts z = 5.7 display 
AB magnitudes z' ~ 24.5 — 25 (Hu et al. 2004). The additional distance 
modulus between a redshift z = 5.7 and z = 8.8 results in a relative 
dimming term of 1 magnitude. 

(2) Optical imaging data must reach a limiting depth typically 1.5 magni- 
tudes fainter than NIR data. Early-type galaxies located at redshifts 
z ~ 2 can generate a spectral discontinuity between optical and NIR 
continuum bands of amplitude D ~ 1.5 mag. (Stanway et al 2004b). 
Failure to identify the continuum break directly could lead to the mis- 
identification of redshifted [OII]3727 emission in such sources as candidate 
z ~ 9 Lya emission. 
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Following these considerations the Hubble Deep Field South (HDFS; Williams 
et al. 2000) Wide field Planetary Camera 2 (WFPC2) apex pointing (a = 
22 /l 32 m 55 s .64, 5 = -60°33'05".01, J2000) was selected as the target field in 
order to exploit the high quality of optical to NIR image data available for the 
field. In particular, the combination of HDFS WFPC2 and Very Large Tele- 
scope (VLT) Infrared Spectrometer And Array Camera (ISAAC; Moorwood 
1997) observations of the field provide images to typical depths AB=28 and 
AB=26 in optical and NIR bandpasses respectively (Labbe et al. 2003). 

VLT/ISAAC further provides a suitable combination of narrow and broad 
band filters for the ZEN survey: the NB119 filter (see Figure 1) is well placed 
within the response of both the J and J s NIR bandpasses. In addition the 
NB119 filter samples a spectral region free from both strong atmospheric emis- 
sion and absorption features. The NB119 filter has an effective width of 89. 5A 
and is centred at a wavelength 1.187/xm (corresponding to the location of Lya 
emission at a redshift z = 8.76). While the NB106 filter samples similarly 
"dark" sky regions, it is not well placed with respect to the broad J-band 
filters - complicating the estimation of NB excess for individual sources. 




1 1.1 1.2 1.3 1 1-1 1-2 1-3 



Fig. 1. Left panel: Filter response functions for the NB106, NB119, J, and J s 
VLT/ISAAC filters. Right panel: The central wavelength locations of the NB106 
and NB119 filters (vertical red dotted line) are compared to the NIR night sky 
emission spectrum (solid back line; Rousselot et al. 2000) and the atmospheric 
transmission measured at two epochs to illustrate temporal transmission variations 
(red and blue solid lines). 
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3 The ZEN1 data set 



Narrow J-band observations of the HDFS WFPC2 pointing were obtained 
during ESO Period 69 (May 19th to September 17th 2002) employing the 
VLT/ISAAC facility. The total NB data set consists of 420 x 300s spatially 
dithered exposures. Imaging data were a) corrected for varying pixel response 
using twilight sky exposures, b) sky-subtracted having masked array regions 
containing objects detected above a specified ADU level, c) corrected for both 
high- and low-frequency spatial artefacts, d) shifted to a common pixel scale 
and coadded using a suitable pixel weighting and rejection algorithm. The 
applied data reduction techniques are broadly similar to those described in 
Labbe et al. (2003) for deep J S HK S observations of the HDFS. A more com- 
plete description of the NB data reduction and analysis can be found in Willis 
and Courbin (2005). 




Fig. 2. Greyscale NB119 image of the HDFS. The image orientation is North up 
and East left. The field geometry of the NIR broad band data set is identical to the 
NB image. The field geometry of the HST WFPC2 field is indicated. Note that the 
Planetary Camera (shaded region) data does not contribute to the final catalogue. 
The total field area contributing to the NB excess catalogue is 4 arcmin 2 . The 
location of the two NB excess sources identified from the catalogue are indicated 
(see text for details). 
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The photometric zero point for the NB119 image was computed employing 
5" diameter aperture photometry of 47 bright, isolated sources common to 
the NB119 and FIRES J s -band images. Assuming that the SEDs of these 
calibration sources display no strong discontinuities, the J s -band flux density 
provides an accurate estimate of the flux density within the NB119 filter, i.e. 
J s — NB = 0. Source detection and photometry was performed on the NB119 
image and the 0."7 diameter detection apertures were transformed to the 
FIRES astrometric system to compute corresponding U 300 B4 50 Vq qIsuJ s HK s 
magnitudes. 

The flux completeness limit of the NB119 and J s images was estimated by in- 
troducing and recovering artificial unresolved sources within each field. Adopt- 
ing the 90% point source recovery threshold as the limiting magnitude in each 
band generates magnitude limits of NB < 25.2 and J s < 26.2 and corre- 
sponds to integrated signal-to-noise ratios (SNR) of 13.8 and 10 respectively. 
The magnitude limit NB = 25.2 corresponds to a total flux integrated across 
the NB filter of F^b = 3.28 x 10~ 18 ergss _1 cm -2 . In common with Labbe et 
al. (2003) we present 0."7 aperture photometry throughout this paper. The 
correction required to convert photometric measures computed in 0."7 aper- 
tures to 5" apertures (which we assume to be 'total' measures) was determined 
to be 0.7 magnitudes via analysis of bright stars in the NB image. 

Figure 3 displays NIR narrow band excess, expressed as J s — NB, versus 
NB magnitude for all sources extracted from the total area covered by the 
optical and NIR data for the HDFS field. Note that all NB detections have a 
corresponding J s detection. Sources displaying J s — NB > 0.3 and NB < 25.2 
are flagged as potential z ~ 9 emitting galaxies. 

Two sources are identified by the above selection criteria, HDFS-352 and 
HDFS-768 - where the identification numbers refer to the HDFS photometric 
catalogue of Labbe et al. (2003). The photometric spectrum formed by the 
UzoqB^qVqqqI^JsHKs plus NB photometry of each object is displayed in 
Figure 4. In each case the addition of deep optical photometry is sufficient 
to exclude each NB excess detection as a potential z ~ 9 source: each NB 
excess object is detected in all optical bands and the photometric redshift 
of each source (Rudnick et al. 2001) - z phot = 1.54 ± 0.06 (HDFS-352) and 
Zphot = 0. 76t(j]2o (HDFS-768) - is consistent with a NB excess arising from 
redshifted H/54861 and Ha 6563 emission respectively. Although no candidate 
z ~ 9 sources are detected within the HDFS survey area, the confirmation of 
two sources whose photometric redshifts are consistent with the observed NB 
excess arising from narrow line emission lends support to the assertion that 
the adopted selection criteria identify faint, narrow emission line galaxies. 
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Fig. 3. Distribution of NB excess J s — NB versus NB magnitude for all objects in 
the NIR field area of the HDFS (solid points). The solid horizontal line indicates 
J s — NB = 0. The vertical dotted line indicates the selection criterion NB < 25.2. 
The horizontal dotted line indicates the selection criterion J s — NB > 0.3. The 
dashed curve indicates the predicted J s — NB uncertainty (3a) as a function of NB 
magnitude returned by the completeness analysis. 
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Fig. 4. Low resolution spectral energy distributions formed by the multi-band pho- 
tometry of the two NB excess objects identified within the HDFS. The clear detec- 
tion of each source within optical passbands precludes each source as a candidate 
z ~ 9 source. 



4 A ZEN paradox: what is the significance of detecting no Lya 
emitters at z = 9? 



The Lya emission line selection function generated by the NB excess search 
technique is displayed in Figure 5 in terms of the co-moving volume sampled 
as a function of Lya emission luminosity for three values of the rest frame 
velocity width of the Lya-emitting source. For the case where the rest frame 
velocity width of putative z ~ 9 sources is a v = 50kms _1 , the NB excess 
survey area samples a co-moving volume of 340 h~ 3 Mpc 3 to a Lya emission 
luminosity of L Lya > 10 43 ftr 2 ergss _1 . 
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Fig. 5. Co- moving volume sampled as a function of Lya emission luminosity. Three 
values of the rest frame Lya velocity width are displayed; a v = 20 (solid line), 50 
(dashed line) and 100 (dotted line) kms -1 . 



The uncertain physical state of the IGM at z ~ 9 complicates a detailed 
comparison of model Lyo-emitting populations to the region of luminosity and 
volume parameter space sampled by the ZEN1 study (e.g. see Santos 2004). 
A more empirical approach is to consider known Lyo-emitting sources and to 
transpose their properties to redshift z ~ 9 . Hu et al. (2004) (hereafter H04) 
present a search for Lyo-emitting galaxies at z ~ 5.7 employing a narrow 
band filter centred on the wavelength 8150A. Though the H04 study does 
not present the highest redshift Lyo sources currently known (c.f. Kodaira 
et al. 2003), their study is based upon near complete spectroscopic follow- 
up observations - 19 out of 23 observed sources (from a total sample of 26 
candidates) were confirmed as Lyo at z — 5.7 - and the overall design of their 
study is very close to the ZEN1 survey. H04 report a surface density of Lyo 
emitting galaxies of 0.03 arcmin~ 2 to a Lya flux limit of 2 x 10~ 17 ergs s _1 
cm -2 f rom a total areal coverage of 700 arcmin 2 . 
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Within our adopted cosmological model, an unresolved source at a redshift 
z = 8.8 is 2.7 times fainter than an identical source viewed z = 5.7. A further 
factor of 0.5 must be applied to account for the different fraction of total flux 
measured by the 3" diameter apertures applied by H04 and the 0."7 diameter 
apertures employed in the current study. The H04 flux limit therefore corre- 
sponds to an approximate limit of 3.5 x 10~ 18 ergs s _1 cm~ 2 when transposed 
to a redshift z = 8.8 and employing 0."7 diameter apertures, i.e., assuming 
that no additional evolution occurs, the z = 5.7 sources observed by H04 are 
sufficiently bright to be observed within ZEN1. Applying the surface density 
of confirmed z = 5.7 Lja emitting galaxies to the areal coverage of the current 
NIR survey, indicates a probability to detect az~9 galaxy of 0.12 - assum- 
ing no evolution between the two redshifts. This 'volume shortfall' indicates 
that the current NIR selected survey will have to be extended by more than 
eight times the currently sampled area in order to realistically probe the no 
evolution scenario. 

The present ZEN1 survey has successfully demonstrated that the sensitivity 
required to detect putative z ~ 9 Lja emitting galaxies emission has been 
achieved for the case where no evolution occurs between redshifts z = 5.7 and 
z = 8.8. However, it is clear that, in order to place stronger limits on the 
space density of such sources, the areal coverage of the current study must be 
extended. 



5 Continuing efforts 

A number of observational programs are underway with the aim of detecting 
z > 7 Lya-emitting galaxies using the narrow J-band technique described 
above (see contributions from Smith et al. and Bland-Hawthorn et al. in these 
proceedings). We are currently undertaking ZEN2, a NB119 survey of three 
massive lensing clusters - AC114, A1869 and A1835. The gravitational lens 
effect provides a typical magnification across each ISAAC field of three and 
our observations have been designed to reach the same unlensed luminosity 
and volume sensitivity as achieved with ZEN1. Of particular interest however, 
are the regions of each cluster field that correspond to critical lines (regions 
of very high magnification) for sources at z ~ 9 . These limited areas in each 
ISAAC image provide exceptionally deep sightlines and offer an important 
'keyhole' through which to detect exceptionally faint emitters. 

Therefore, with the completion of ZEN2 and existing NIR NB surveys (most 
notably DAZLE; Bland-Hawthorn et al. in these proceedings) the pursuit of 
z > 7 Lya-emitting galaxies will have reached the stage where simple and 
robust galaxy evolution models can be tested directly. 
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